one-dimensional Ising model built on small-world networks: competing dynamics 
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In this paper, we offer a competing dynamic analysis of the one- dimensional Ising model built 
on the small- world network (S WN) . Adding-type SWNs are investigated in detail using a simplified 
Hamiltonian of mean-field nature, and the result of rewiring-type is given because of the similarities 
of these two typical networks. We study the dynamical processes with competing Glauber mechanism 
and Kawasaki mechanism. The Glauber-type single-spin transition mechanism with probability 
p simulates the contact of the system with a heat bath and the Kawasaki-type dynamics with 
probability 1— p simulates an external energy flux. By studying the phase diagram obtained in the 
present work, we can realize some dynamical properties influenced by the small- world effect. 
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Since Watts and Strogatz proposed the small- world networks (SWN) [1], which are believed to catch the essence of 
many network systems in nature and society, a large number of further works have appeared (see [2-5] for review). 
Researchers are interested in investigating the properties of various models and processes on SWNs. Recently, Zhu 
and Zhu. successfully introduced the SWN effect to the critical dynamics of the spin system [7] , and thus extended 
the investigation to the dynamic properties of spin models. In the 1960s, the dynamic behavior of the Ising model 
was successfully described with the Glauber [8] and Kawasaki [9] mechanisms. From then on, great progress has been 
achieved based on these two mechanisms. Besides, an interesting problem has been attracting much attention, i.e., 
the competing Glauber-type and Kawasaki-type dynamics [10-12]. In the competing mechanism, the Glauber-type 
3 ■ dynamics is given probability p, while the Kawasaki- type one has probability 1 — p. Zhu et a/.'s recent work [13] has 
studied small-world network effect in the competing dynamics on the Gaussian model. Some meaningful results has 
been obtained, but due to the requirement of the convergence of the integration, they were not able to get the full 
phase diagram. In this paper, we investigate the competing dynamics of the Ising model considering the small-world 
network effect, and we obtain the full phase diagram and the competing dynamic behavior. By this work, we can 
further understand the influence of the SWN effect and highlight the disparities between the dynamic mechanisms. 

In the construction of SWN with a certain probability of introducing long-range links, we will end up having a 
whole set of possible realizations. Thus the theoretically correct way of treating dynamic systems built on SWN should 
\ involve three steps: First we have to make a full list of all the possible realizations and point out the probability 
of each one of them. Second, we treat the problem independently on each network. Finally, we give the expected 
value with all these results. Although being conceptually straightforward, it is very cumbersome even for the simplest 
one-dimensional Ising model. It has been suggested that the spin system on the SWN as a whole has mean-field-like 
effect on individual spins due to the long-range links [6,7]. Naturally, a simplified method of mean-field nature was 
presented by Zhu et al. [7] to describe the kinetic spin system built on SWN. According to this simplified method, all 
possible networks are deemed as a single one. Then, the effective Hamiltonian of a spin-lattice model built on such a 
network is defined as the expected value over all possible realizations, 
i-^j ■ In the present work we study two specific networks: In a one-dimensional loop, for example, (1) each randomly 
selected pair of vertices are additionally connected with probability p A ; and (2) the vertices are visited one after 
another, and its link in the clockwise sense is left in place with probability 1 — p R and is reconnected to a randomly 
selected other node with probability p R . Networks of higher dimensions can be built similarly. We shall refer to the 
J> , first model as adding-type small- world network (A-SWN) and the second one as rewiring-type network (R-SWN). 

We discuss the problem using the simplified method, and give the effective Hamiltonian first. For the one- 
dimensional (ID) Ising model built on A-SWNs and R-SWNs, the effective Hamiltonian can be written as, respectively 
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where (3 — l/k B T and K = (3 J. The case of K > 0, (J > 0), corresponds to the ferromagnetic system. 

Various dynamic processes in critical phenomena are believed to be governed by two basic mechanisms, i.e., the 
Glauber-type with order parameter nonconserved and the Kawasaki-type with order parameter conserved. Their 
combination, namely the competing dynamics, gives 



dt 



P({a},t)=pG me + (l-p)K TI 



(3) 



d_ 

dt 



q k {t) = V Q G k + (1 - p)Q* qi (t) = J2 °iP (M, *) • 

M 



(4) 



This dynamic competition has also been receiving attention [12,13]. Here, pG me (or pQ'j?) denotes the Glauber-type 
mechanism with probability p and (1 — p) K me [or (1 — p)Q k \ denotes the Kawasaki-type mechanism with probability 
1—p; they are determined, respectively, by the Glauber-type single-spin transition probability Wi (<7j — > a{) [14] and 
by the Kawasaki-type pair-spin redistribution probability Wji(ajai — > <7j<7;) [15]. 

In their original form, the Glauber-type dynamics and the Kawasaki-type dynamics both favor a lower energy state. 
However, the competing dynamics is usually used to describe a system in contact with a heat bath while exposed to 
an external energy flux. Naturally one requires a competition between one process favoring lower system energy and 
the other one favoring higher system energy. Usually, the Glauber-type mechanism is used to simulate the contact of 
the system with a heat bath and it favors a lower energy state. On the other hand the Kawasaki-type mechanism can 
be modified to simulate an external energy flux that drives the system towards a higher energy state. This can be 
achieved by switching (3 to — (3, or K — (3 J — J/KbT to — K, and modifying the redistribution probability accordingly. 
This means that the competition between the Glauber-type mechanism and the Kawasaki-type mechanism is actually 
a competition between ferromagnctism and antiferromagnetism. Probing the competing behavior of the ID Ising on 
A-SWNs and R-SWNs is certainly a problem of interest. The equation of evolution of the local magnetization is given 

by 
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The first and the second term correspond, respectively, to the part of the Glauber's dynamics with probability p, and 
the part of the modified Kawasaki's dynamics (only switching K = (3 J to —K) with probability I— p. Our calculation 
will focus on A-SWN, while the result of R-SWN will be given straightforwardly. 
For the Glauber-type dynamics 
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and for the Kawasaki-type dynamics 
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where the expressions of Wi (<Tj — > &i) and Wji(ajai — ► OjGi) can be found in Rcf. [7] (but differently from that paper, 
we should switch (3 to — (3 for Wji). Three important combined terms in Eqs. (6) and (7) are calculated, for the ID 
Ising model, to be 
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Now we turn to determine the system behavior by studying the tendency of evolution of the following order 
parameters: 
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Obviously, a state with both vanishing M(t) and M'(i) corresponds to the disordered paramagnetic phase; a state 
with nonvanishing M(t) and vanishing M'(t) corresponds to the ferromagnetic phase; and, in an antifcrromagnetic 
phase we will find vanishing M{t) and nonvanishing M'(t). If both order parameters are nonvanishing, this phase 
cannot be simply identified as a ferromagnetic or paramagnetic phase, but can be tentatively named as a heterophase. 

The cases of A-SWN will be discussed in detail in the following, and we will determine the tendency of system 
evolution under a small perturbation M(t) — ► 0. We then give the results of the case of R-SWN. 

For the ID Ising model built on A-SWNs, the effective Hamiltonian of the system is given by Eq. (1). From Eqs. 
(6)- (8), we can easily obtain 
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The tendency of evolution of the order parameters M(t) and M'(t) can be demonstrated by Eqs. (9) and (10) when 
the system undergoes a small perturbation. 

(1) When p A = 0, it means that no long-range link exists, and 



M(t) = M exp 



-p i 



tanh 2K \ 
tanh 2^0 J 



(11) 



M'(t) = M^exp 



P U + 



tanh 2K 
tanh 2K { }. 



+ (l-p) 2 



tanh 2K 
tanh 2K ( > 



(12) 



where tanh2_ft'" = 1 corresponds to the critical point of the one-dimensional (ID) Ising model without the SWN 
effect (K° -> oo, T c ° = 0). When K < K°, (T > T c ° = 0), Eqs. (11) and (12) show that M(t) and M'(t) are both 
approaching zero exponentially due to the fact that tanh2if < 1, and thus we can reach the conclusion that, by 
whatever amount one increases the energy flux, the system will stay in the paramagnetic phase at arbitrary finite 
temperature. When K -> K°, M'(t) ->■ 0, and 
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the critical slowing down of the order parameter M(t) will appear at the critical point K%. 

(2) Now a small portion of adding-type long links are introduced to the system, the system behavior in this case 
can be described by Eqs. (9) and (10). Obviously, M'(t) is approaching zero exponentially at any temperature and 
any p A . From Eqs. (10) and (12), we can see clearly that the rate of M'(t) approaching zero is faster than that in 
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the regular network. Different from M' (t), the evolution tendency of the order parameter M (t) depends on both K 
and p A . We can obtained the critical point by the following equation: 



tanh 2K A + (N-1) p A K A (l - ^ tanh 2 2K A ^j = 1. 



To give an example, if we suppose p A = 1/N, we can get t&nh.2K A \ pA ^ 1 ^ N — 0.6809, or K A \ pA ^ l / N = 0.154. Relative 

to K A , when t — > oo, M(t) — > for K < K A , M(t) ^ for K > K A , and M(t) experiences critical slowing down for 
K — > i'Q 4 . So, combining M (t) and M'(t) we can conclude that: (2. a) For the case of K < K A , the system stays in 
the paramagnetic phase; (2.b) For the case of K — > if^, the system shows the critical slowing down; and (2.c) For 
the case of K > K A , the system stays in the ferromagnetic phase. 

However, when p = 0, we cannot identify the system simply as ferromagnetic or paramagnetic. Because in this 
case, it depends on the initial state. If Mo ^ 0, the system will stay in ferromagnetic, otherwise the system will be 
paramagnetic. 

The phase diagram is shown in Fig. 1(a). 

For the ID Ising model built on R-SWNs, the effective Hamiltonian of the system is given by Eq. (2). With 
analogous calculation, one can get the equation of the critical point 

tanh [2K R (1 - p R )] + 2K R p R jl - \ tanh 2 [2K R (l - p R )] } = 1, (13) 

and the time-evolution of the orders parameters M(t) and M'{t) 
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(1) When no rewiring long link exists, i.e., p R = 0, the evolution is the same as the case of p A = 0. 

(2) When a small portion of rewiring- type long links is introduced to the system, M'(t) is approaching zero 
exponentially at any temperature and any p R . The decay rate of M'{t) is also faster then before. Different from M' (t), 
the evolution tendency of the order parameter M (t) depends on both K and p R . The calculation of the critical point 
is similar to the A-SWN one. When p R = 0.1, we can get the critical point, tanh [2K R (l - p R )] \ pR=Q 1 = 0.90222, 

or K*\ pR=0 1 = 0.82446. Relative to K?, when t -> oo, M(t) for K < K?, M(t) ^ for K > K*, and M(t) 

shows critical slowing down for K — ► K R . Combining M (t) and M'{t) we can conclude that: (2. a) For the case of 
K < K^, the system stays in the paramagnetic phase; (2.b) For the case of K — > K^, the system shows the critical 
slowing down; and (2.c) For the case of K > K^, the system stays in the ferromagnetic phase. 
The phase diagram is shown in Fig. 1(b). 

In this paper, we analytically study the dynamic properties of the ID Ising model built on small-world networks. 
Two typical SWNs are investigated, the adding type and rewiring type. As is generally known, the ID Ising model 
on the regular lattice does not show continuous phase transition at any nonzero temperature. However, if the SWN 
effect is introduced, critical phenomena appear in the ID Ising model. With competing dynamics, as long as p ^= 0, 
the phase diagrams are separated into two regions. Below the critical temperature, the system will get into the 
ferromagnetic phase, while above the temperature, the system will get into the paramagnetic phase. The critical 
temperature is independent of the competing probability p. Different from the 2D Ising model [10] and the Gaussian 
model [13], the ID Ising model built on SWNs does not show antiferromagentic phase at any temperature and any 
competing probability p. 

As we have seen above, the Ising model shows critical phenomena on both A-SWNs and R-SWNs. This is because 
random links introduce long-range interactions. It is reasonable that the system will exhibit long-range order at finite 
temperature. Furthermore, the more extra links, the higher the critical temperature. Thus on R-SWNs the Ising 
model has a maximum critical temperature for there is a maximum number of the random links (it is also expected 
that the simplified method with an effective Hamiltonian will fail when most regular links are rewired). However, 
long-range interactions do not lead to antiferromagnetic order at any competition probability, but instead, the long 
range links make any antiferromagnetic order decay faster. This is because the long-range interaction here is random, 
while antiferromagnet needs ordered long interactions. 
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Caption of figures 



Fig.l. The phase diagrams of II? fsing model built on SWN with competing dynamics, Glauber-type with proba- 
bility p, and Kawasaki-type with probability (1 — p). (a) On the A-SWN structure: p A ~ 1/N; (b) On the R-SWN 
structure: p H = 0.1. In which, tanh2iir° = 1 corresponding to the critical point of the ID Ising model without the 
SWN effect (K° — ► oo, T c ° = 0). 
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(a) A-SWN: p A ~l/N 



(b) R-SWN: p R =0.1 
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